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A NONUNEAR THEORY FOR SONIC -BOOM CALCULATIONS 


IN A STRATIFIED ATMOSPHERE 

By Michael Schorling* 

Langley Research Center 

SUMMARY 

I 

The exact solutions to the equations of gas dynamics are given with respect to the 
axis of slender lifting bodies in a stratified atmosphere. The boundary condition is sat- 
isfied by using slender-body theory. The solution predicts the magnitude of the pressure 
rise of the sonic boom and estimates the nonlinear effects in the vicinity of the cutoff 
point. 


INTRODUCTION 

The theories of Friedman, Kane, and Sigalla (ref. 1) and Hayes, Haefeli, and 
Kulsrud (ref. 2) on sonic-boom propagation in a stratified atmosphere have the advantage 
of being able to treat nonsteady flight maneuvers. However, they have the disadvantage of 
not being fully nonlinear theories. Nonlinear effects such as focusing at the caustic are 
now of great interest. Here a theory is presented in which the exact nonlinear equations 
for supersonic flow are solved for large distances from lifting bodies. This paper is an 
extension of the author's previous papers (refs. 3 and 4) to a stratified atmosphere in 
which the speed of sound changes with the altitude of flight. The pressure p and the 
density p are assumed to obey the hydrostatic law dp = -pg dz. The gas is considered 
to be nonviscous, steady, homentropic, and homoenergetic. The shock front should be 
attached at the nose of the slender pointed body. 

The theory presented can be used to calculate the gas properties in the far field of 
nonaxisymmetric bodies by employing area-rule concepts. 

SYMBOLS 

A cross-sectional area of body 

A' first derivative of cross-sectional area of body 

♦ 

NRC-NASA Resident Research Associate. 




speed of sound 


a(r) 


speed of sound at distance r from the body 


a(z) 


speed of sound at altitude z 


specific heat at constant pressure 


v specific heat at constant volume 

3x 8 8xd 


D '-'./V ^ '-'SV v r\ rv 

r X t 8 r - X r 8 1 

84 3r 8r 84 ? r ^ 


Dt = — = 8t 
? 84 $ 


_ 8x __8_ 8x_ _8_ 

“ 04 dtp 8 tp 04 


,1, - — — - — — = X^8^ - X^8^ 


E 


= 2£ 


e 1 = grad x 1 


e" 2 = grad x2 




/ Cartesian base vectors 


e 2 = grad x^ 


p(j) 

m 




g 

g 


function of integration 
function of integration 
acceleration due to gravity 
volume force due to gravity 


H (1) (A)I 

) functions dependent on the cross section (see eqs. (44)) 

H (2) (A) 


h specific enthalpy 


2 


Mach number 


generatrix of Mach cone 

number of degrees of freedom of the gas (air: n = 5) 

generatrix of the wave normal cone 

pressure 

distance from the body at which the boundary condition is satisfied 

radial distance from the body 

entropy 

locally dependent base vectors 
time 

gravitational potential 
magnitude of the velocity 
velocity vector 

components of the velocity vector in a locally dependent basis (i = 1, 2, 3) 
abscissa of the characteristic surface 
abscissa of the shock front 

xp > Cartesian coordinates 

■J \ 


altitude of flight 


a 


angle of attack 


(3 Prandtl factor, \j - 1 

y ratio of specific heats 

A = w<°) 2 - a<°>(r) 

Ap = P - Poo 

r) Prandtl transformation, /3r 

3 angle of inclination of the streamlines 


A = 




( 0 )' 


a<<) 


1/4 


a«»(r) 


n 


2\/r 


li Mach angle 

£ characteristic variable 

^ = (£,r,i//) independent variables in locally dependent basis (i = 1, 2, 3) 


n = 


[a 


a^(z)J 


-.(n+l)/4 


w 


\2 / n \ 2 "1 


1/4 


(or _ a (o) 


(z)j 


density 


S = a(°)(r) 


ri = (w ,d ,(p+ip) where i = 1, 2, 3 


<P 


azimuth angle 


4 


<p 


potential function 


\p cylindrical variable (azimuth angle) 

co variable of integration 

Subscripts: 


00 state of undisturbed flow 

gr values at the ground 

Harmonic (or Fourier) numbers are denoted by numbers in the subscript position or 
by the subscript m. The order of magnitude is indicated by numbers in parentheses in 
the superscript position or by the superscript (j). The Greek letters A, a, and v 
used as subscripts and superscripts denote the Einstein summation convention. A comma 
preceding a subscript denotes differentiation with respect to that subscript. 

ANALYSIS 


Basic Equations 

The gas considered here is assumed to be in equilibrium. There are no heat or 
mass sources. To describe the state of the gas, it is sufficient to calculate the vector of 
velocity, the density, and the pressure. This can be done by means of the equations of 
continuity, momentum, and energy. From these equations it may easily be seen that 

f = ° <« 


As the entropy is assumed to be constant in the overall undisturbed field and, according to 
equation (1), it does not change along streamlines, the entropy is constant in the disturbed 
field as well. That is, the gas is homentropic. ■ 

Now the change of energy is to be considered. The equation of momentum is 


(w • grad)w = - -i grad p + g 

r 

Combining equation (2) with the identity 

— — w 2 _ 
(w . grad)w = grad — — w x curl w 

Li 


(2) 


(3) 


5 



gives 


w 2 _ i _ 

grad — - w x curl w = -- grad p + g 
2 P 


where 


( 4 ) 


= w 


w 


As the volume forces due to gravity are irrotational, 
g = grad U = (0,0, -g) 

Projecting equation (4) onto the direction of the velocity vector w yields 



as w x curl w = 0 along the streamlines (Beltrami flow is not considered here). Equa- 
tion (5) can be integrated for a given relation between the pressure p and the density p. 
According to thermodynamics, 


P 


P(p,s) 



( 6 ) 


where the bar over symbols indicates an arbitrary initial state of the gas. For homen- 
tropic flow, equation (6) is rewritten as p = p y . Const. Thus, the integration of equa- 
tion (5) yields 

+ gz + h = H (7) 

Ct 

where H is the total enthalpy. Along the streamline the total enthalpy is a constant. 
Since the undisturbed flow field was assumed to be homoenergetic, the disturbed flow field 
is homoenergetic as well, according to equation (7). As long as no shocks are taken into 
account, the result of these considerations is that the overall flow field is homoenergetic 
and homentropic. Far from the body weak shocks are expected. For weak shocks the 
entropy increases as the third power of the pressure coefficient. Thus, so long as only 
weak shocks are considered, the flow can be assumed to be homentropic. 
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The basic equations can now be formulated. According to equation (1) the flow is 
irrotational. This condition is expressed by the equation of irrotational flow: 

curl w = 0 (8) 


The equation of continuity is 


— + p div w = 0 
dt 


>1 


0 ) 


By use of the relations p = p(p) 


and 



equation (9) can be rewritten 


1 dp 
P dt 


+ a 


2 div w = 0 


For steady flow this equation can be expressed in terms of the speed of sound, the veloc- 
ity, and a term due to gravity: 

a 2 div w - ww • grad w + w • g =0 (10) 


The last basic equation is derived by considering the position of the characteristic 
surfaces. Assume that the quantity | denotes a characteristic variable along the gen- 
eratrix n of the wave normal cone. The unit vector in along the generatrix of the 

Mach cone is perpendicular to n by definition; that is, in • rf = 0. As n = ■ - - 

I grad || 

or, according to the definition of the speed of sound, (w . n ) = -a, the last basic equation 
is written in terms of the speed of sound, the velocity w, and the characteristic variable 


9 “*-2 

2 r grad | • grad | = (w • grad |) (11) 

Change to a Characteristic Coordinate System 

The basic equations have been written in an invariant form. Now a coordinate sys- 
tem must be selected to carry out the solution. In a Cartesian frame, a point has the 
cylindrical coordinates xl = x, x 2 = r cos \j/, and x 2 = r sin \ f/. However, here a ^locally 
dependent basis will be used, in which the velocity vector w has to be expressed. 
According to figure 1, the velocity vector w is given by 
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w = w 


cos 9 e"* + sin 9(cos <p s’ + sin <p t) 


(12) 


where 


s = cos 


if/ e 2 


, —3 
sin i// e 


— _o —3 

t = -sin ip e + cos ip e 

The velocity vector is defined by its magnitude w, by the angle of inclination 9 of the 
streamline, and the azimuth angle cp + ip. If w, 9, and cp + ip are considered to be 
new variables of the velocity vector, new components of the velocity are defined as 

w i= ^ ri = (w, 9, <?+<//) (i = 1, 2, 3) (13) 

&T 1 


Up to now the unknown functions have been expressed by the cylindrical coordinates 
w = w(£,r,ip), i; = £(x,r,ip), and a = a.(x,r,ip). Here the problem is to calculate the posi- 
tion of the characteristic surface rather than to determine the characteristic surface 
itself. So, following the Poincare- Lighthill-Kuo (PLK) method, the dependent variable £ 
is exchanged for the independent variable x. The new set of independent variables is now 
£, r, and ip. The velocity vector w, the speed of sound a, and the abscissa x are to 
be determined as functions of these variables. That is, defining ^ = (£,r,i//) where 
i = 1, 2, 3, determine w = w(^), x = x(£j), and a = a(£j). The frame of reference of 
these coordinates must be constructed and the differential operators determined in these 
given directions. The final results are (see ref. 3), 


x^ grad = e X D| + sD r + ^ ?D^ 

Q — ► 

curl w = grad t x w CT 
div w = grad r a . w CT 


(14) 


where D^, D r , and are abbreviations for the operators defined by the equations 

= 9^, D r = x^9 r - x r 9^, and D^, = x^9^ - x^9|. Equations (12) to (14) have to be 
introduced into the basic set of differential equations (eqs. (8), (10), and (11)) to carry out 
the solution in the given frame of reference fe^s", t). 
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— —2 

Since g = -ge , the equations are written as functions of the coordinates 

«i = (4,r,*): 

sin B sin <p D r w + -p sin 3 cos cp D^w - w cos B sin <p D r 3 + ^ cos B cos cp D^B 
- w sin a cos <p D r ^ - ^ sin B sin cp D ^cp - ^ sin a sin cp x^ = 0 ' (15a) 


sin B sin cp D^w - -p cos B D^w + w cos ,9 sin cp ^ sin B D ^B + w sin 3 cos cp D^c p = 0 

(15b) 


-sin 5 cos cp D^w + cos 3 D f w - w cos 3 cos cp D ^B - w sin 5 D r 3 + w sin B sin cp D^cp = 0 


^a 2 - w 2 ) 


1 


cos B D^w + sin B I cos cp D r w + - sin B D^w, 


o 

+ a^w 


-sin B D gB 


1 


+ cos B |cos cp D rl 9 + - sin B D ^B 
- x^gw sin 5 cos ip = 0 


+ a 2 w sin B 


-sin cp D r <p + - cos cp D ^cp 


(15c) 


(15d) 


w 2 fcos 5 - sin B cos cp x r - sin B sin cp i x^j = a 2 (l + x r 2 + x^ 2 


w 


2 + na 2 + 2gz = Constant 


(15e) 

(15f) 


The last equation - the energy equation — gives the interrelation between the speed of 
sound and the velocity. Although this set consists of six equations for the five unknown 
functions w, B, cp, a, and £, all of them must be used. The equations curl w = 0 
are dependent on each other and have, together with the equation of continuity,, the trivial 
solution w = grad <p, since curl • grad cp = 0. 

i 

Solution by a Perturbation Approach j 

1 

The system of nonlinear partial differential equations (15) is solved for large dis- 
tances by a perturbation method. If the unknown functions are assumed to have the form 
of a perturbation series 


9 


+ • • 


(16) 


c = C^ + + C^ 

that is valid for large distances, the nonlinear differential equations can be split up into 
orders of magnitude. The equations for the different orders of magnitude can be easily- 
solved in a step-by-step procedure. The series (eq. (16)) does not make any statement 
about the values of the orders of magnitude. It only expresses the assumption that for 
large distances r the zeroth-order term is larger than the first-order term, and so 
forth. The special properties of this series will be determined in this paper. There is 
some a priori knowledge about the zeroth order which is put into the calculation. For 
large distances r, to zeroth order 

w^°) = Constant a^ = Constant 3^ = 0 = Constant x^ = 0 (17) 

Contrary to the previous paper (ref. 3), x^ is now dependent on r. This complicates 
the calculations greatly. Furthermore, for large distances the following assumptions are 
made concerning the derivatives: 

8 r 9 | ~ 9 ^ 9 r « a £ 9 r « 8 ^/ ( 18 ) 


but 



According to the linearized two-dimensional equation of characteristics x = £ + /3r, x^ 
has been set equal to 1. Equations (16) to (18) must be introduced into equations (15) to 
separate the equations for the different orders of magnitude. 

The final differential equations are, to zeroth order, 

w ( o ) 2 _ a (0) 2 

and to the first order, 


1 + x 


(or 

r J 


= 0 


(19a) 




= o 


(19b) 


w<°Vl> + xV = 0 

£ r £ 


(19c) 
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x<°>w£> + w(°^ = 0 (19d) 

a (0) 2 x (0) x (l) . .,(l) x (0) w (0) 2 + w (0) w (l) . a (0) a (i)[ 1 ■ x (o) 2 ] (19e) 

w (°) ^V 1 ) ^ = 1 w^ (19f) 

w<°>[> + 1 S (lj] . A(>) . s (l) = 0 (19g) 

L a ® 2 

na^a^ + w^V 1 ) = 0 (19h) 

While w^) is a constant in the whole field, a^ is dependent on the altitude of 
flight z. By means of the integrated hydrostatic law, 

na 2 (z) + 2gz = na 2 (20) 


the speed of sound can be determined for an arbitrary height or given distance from the 
body. In equation (20), z is the altitude of flight, a(z) is the speed of sound at the alti- 
tude z, a is the speed of sound at the ground, and n is the number of degrees of 
freedom of the gas. The local speed of sound a(r) at the distance r away from the 
body is, according to figure 2, 

a 2 ( r ) = a gr “ ( z + r cos (21) 

To indicate the zeroth-order form of equation (21), this equation is rewritten as 
2 2 

a^(r) = a 2 r - E(z + r cos xj/) = a^(z) - E cos \]s • r (22) 

where E denotes a factor depending on gravity, E = 2g/n. 

By developing a Fourier series in xj/ the number of independent variables in equa- 
tions (19) can be reduced from three to two. The proposed Fourier series is 
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(23) 


C (j) (£,r,<//) = £ cjj»>(*,r)e iX * (j g 1) 


A=-°° 


for the variables x, w, 3, and a and 

OO ' ' r • ! 

?>^(l»r,^) = Y, ic Pxh^,r)e lX ^ . 

A=-°° 


(j — ;1) 


(24) 


for the azimuth angle cp. Naturally, the number of equations is increased by this pro- 
cedure from seven equations in three variables to an infinite number of equations in two 
variables. However, for a lifting body of revolution the first two terms of the Fourier 
series will suffice. 

"r. 1 

After having introduced the Fourier series (eqs. (23) and (24)) into the set of equa- 
tions (19), the differential equations are rewritten as 




(0) 2 w<°> 2 -a<°> 2 

‘ r " JO ) 2 


(25a) 


(25b) 


w (o) a a) (o )w a) o 

m,£ r m,| 


(0) a (!) , w(Q) 9 (1) ; C (0) W (1) gW W CQS xp 3 (1) 0 


(0), 


w 


m,r r m r m,r 


,( 0 ) 


2 m 


(25c) 


(25d) 


A<» 

+ w^S^ 

r 

m 

m 

w<°>' 

— 

-i 

i 

M) 3^ 

Vm-v 3 v 


V 

_ 


r m 


j 


(0) 2 J0),(1) . . J1) V (0),„(0) 2 ^ „,(0)„,(1) _ J0)J1)[ + x (0) 2 ' 

r _ 


a x r x m,r “ 3 m x r w + w w m a a m 


(25e) 

(25f) 

(25g) 
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(25h) 


na 


(°) a (l) + w (0) w (l) = o 


m 


m 


Here the subscript m indicates Fourier numbers and the superscripts indicate the order 
of magnitude. The summation index v can be restricted to 0 = k = 2 as the body. . . . 
should be symmetrical with respect to the plane ip = 0 and the first Fourier coefficient 
is assumed to be much greater than the next one, and so forth. 

The solution of equations (25a) and (25b) can be written as - ' 


<P 


( 0 )_ 


Constant 


(26) 


and 


x 


( 0 )_. 1 


E cos \p 




. a»(r) 


w 


(0)‘ 


E cos ip 


|tan‘ 


a<<) 


w 




- tan 


- 1 . 


aKr) 


w 


(or 


a«»(r) 


+ 4 


(27a) 


For E = 0 or cos \p = 0, l'Hospital's rule gives the limit: 


x<°) = r 


w 


(°> 2 - a (0)? r) , 


a (0) (r) 


(27b) 


The integration of equation (25c) with respect to 4 gives 


4 0) w« + w(») S W=w(°)FW(r) 


(28) 


If the integration function F^(r) is set equal to zero, equation (28) becomes the same 
as equation (25e). 

Differentiating this equation with respect to r yields 


i(0)w(l) + x (0) w (l) + w (0) a (l) = 0 
r m. i r . m,r m,r 


(29) 
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Adding equations (29) and (25d) eliminates the quantity w 


( 1 ) . 
m.r' 


23 ^ + 
m,r m 


1 g cos \p x r 
a«»(r) x r 


( 0 ) 


(0) 


= 0 


( 30 ) 


The integration gives 


^(i) = G^\4)r" 1//2 n 
m m vsy 


A 1 / 4 

s (n+l)/4 


(31) 


Here G^(£) is a function of integration which is to be determined by the boundary con- 
dition. For brevity the abbreviations II, A, and 2 are used henceforth: 


n = 


_a«»(z) 


(n+l)/4 




w«» 2 - a<<) 


1/4 


A = w<°> 2 - a«»(r) 


2 = a^(r) 


From equation (28) the first-order velocity component is calculated as 


(!) _ ...(0) G m (g) n 


w m =w 


\fr A l/4 s (n-l)/4 


(32) 


Equation (25g) can be rewritten, by using equations (25h) and (28), as 

-s 2 


(1) _ n + 1 
m,r n 


1 + x 


r (0)‘ 


( 0 )\ 

r J 3 (1) 
% 
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or 


x 


( 1 ) 

m,r 


n + 1 
n 


G in ) (a w(°) 4 n 
'/r A 3/4 2 (n+5)/4 


(33) 


This equation is solved by developing the terms and in a power 

series. 

After the development in a series in terms of r, the integration can easily be car- 
ried out: 

/ 

x4 


X (D = - n + 1 G (1) (t) 

x m n m ™ 


w 


(0) 


a<°>(z) 


v,<°) 2 - a<°)( 2 ) 


2x /F - E cos £ r 3 / 2 
V 6 


■*\ 


Lw«» 2 - a<°>(z) 


n + 5 




+ (E cos xpf r 5/2 
80 


21 


n^ + 14n + 45 


r 


6(n + 5) 


a (°)(z)[w<°) 2 - a<°)(z3 


Lw(»> 2 - a<°>(J 

\ 


/ 


a<°)i 


Ia"' ) (z) 


(34) 


Here the function of integration has been chosen to be zero. The calculations still con- 
tain the quantity n, the number of degrees of freedom of the gas. For air, which obeys 
the model of a molecule with two atoms, n = 5. 

By a development analogous to that of the first-order set of equations, a set of 
second-order equations is written as 


JlUD 

x (0) w (2) + W (0U2) = (1) . "I W 
r ? 1 r ,<°)„,(0) 


Xy 'W' 



(35a) 
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-x!°>wi 2 > + w<°> fl < 2 > + 1 w (0) a (2) . g COS »y») a (2) = + JL±1 |i x <0) 


Aor 




E, 

[ 1 + x r 


UU1)r ''' 4l+ w ( ^(4 0)2 ^E-x( o)2 | 

E 




■ 


(35b) 


x (0) w (2) + w (0) s (2) = w (0)| (l) s (ljl + w (0) (l) a ( 1); . *dU» < 

r ¥ / r LJ — » ,r ' -T ; . ,W (o) 

x r 


,( 0 ) 


1 .1 

•+ — 


x (°) 

v?r 


2 n 


1 + x 


(0)' 


,(0) 


2x< 0 >x< 2 >a< 0 > = 2w(°) v 



r 


2 1 

[uxf^Vl.wW 2 * 

1 

In + l\ 2 

l + x<°) 2 

r 

k )2 -'].4 o,2 f 

\ » / 

x<°> 



L r J 

J 


(35c) 


w(D 2 + 2w<°>w< 2 > + nad) 2 + 2„a<°>a< 2 > = 0 


(35e) 


The formal solution for this set of equations is obtained in a manner similar to that 
for the first-order equations. By investigating the solution for the second-order angle of 
inclination it can be shown that this quantity is really much smaller than the first- 

order angle of inclination. The solutions will be provided to enable an investigation of the 
behavior of the second-order quantities at the caustic, that is, the condition w2 -a2(r)«l 
(However, the second-order calculations have not been included in the sample calcula- 
tions.) The differential equation for the second-order angle of inclination is 


o a (2) + ,(2) 


(o) 2 

g cos \p w v 'E cos 4/ 

_n 2 n (1) G (1) j 

i A E cos 

CM 

% 
i 

2A2 

a ° (r) J 

. 1 

1 r 2 2 7/ 2 

2A 2 2 5//Z 


5A 


22 


9/2 


+ n 2 G^' G^ E cos £ 


m-y v nr 


-11 w 


(0) 2 w (0) 2 


2 _9/2 9 5/2 

? 2A 2 E 


n r 1 ; d« { (- f + n-.2j .'ill— + E^7 3 

Jn m . \ 4 / 5/2 . 1/4*, _3/2 L . 1/ 


r 5/2 A l/4 s r 3/2 l 4A l/4 E 2 , 2 a 5/4 


(E cos |//) 2 

-3 

r l/2 

2A 1 / 4 2 3 


7w 


(0)^ 


21 


i(°) 

F^\r) - F^{r)-£— 
m ' ' m ' ’ (o) 


(36) 
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The left-hand side of this differential equation has the same structure as the analogous 
equation in the first-order calculation (eq. (30)). Hence, the second-order angle of incli- 
nation can be determined as 


^ = nG m ^ i / G J ^ f 

m m-v v j c 


,3/4 


00 r 3/2 S 2 


dr + 2G^A + 




S 1 / 2 

A^r 2 


dr 


+ nAG^^G^E cos ^ 


r / w^ 


5A 


3/4 


0O '2r 1// ' 2 A 9//4 S 2r 1//2 S^' 


dr 


HAG^ flW ^ C0S ^ W f 

m-f v n J c 


<0 »>/ -n + J — _l dr 


^r^A^S 3 2r 1//2 A 9//4 £ 


+ AE cos 


, 1-/2 


00 WrA 1 / 2 ^ 1 / 2 2A 3/ \ 


dr 


+ a(e cos 


r r I -3 ‘ 7w^ 

+ 


2 


21S 


l/2\ 


0 m ^Ua^s 3 / 2 sa 5 / 2 ^ 1 / 2 I 6 A 5 / 2 , 


dr 


(37) 


where 


'.‘ t *. 1 '.5 •<’ 


"■ A 1 / 4 

a. = — 

2S 3 / 2 r 1 / 2 


(o) 

and G v ’ is the second-order integration function. In calculating the integration func- 
tions G^ and G^ , it will be shown that G^ is really much smaller than G^. 
The integration functions are assumed to obey the relation (eq. (23)): 


G (i) = G^(|,^) G^U)e lX ^ 
X 


For completeness, the second-order differential equation for the abscissa of the Mach 
cone will be given- as well. According to equation (35d), 
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(38) 


If the influence due to gravity is set equal to zero, all the previous equations can be 
reduced to the corresponding equations in reference 3. 

As this theory is based on the hyperbolic differential equation of wave propagation, 
the solution to the problem at the cutoff point w = a(r), which is a parabolic problem, 
cannot be obtained. Here, the Mach cone and the wave normal cone degenerate to a plane 
and a straight line, respectively, when the Mach angle p. = 90° (fig. 3). 

Coordinates of the Shock Front 

The differential equation for the coordinates of the shock front in a stratified atmo- 
sphere does not differ essentially from equation (28) of the previous paper (ref. 3). For 
a stratified atmosphere, this equation is written as 


r i 

2 / \ / .2 
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This equation must be developed in orders of magnitude and in Fourier series in in 

the same manner as the set of equations (15). The coordinate x of the shock front is 

s 

calculated to be 


x 


s 




( 1 ) 

m 


(40) 


In this paper the bow shock is the only shock that is considered and calculated. 
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Satisfying the Boundary Condition 

The integration of the system of equations (25) and (35) introduced the functions of 

integration G^(£), which have to be determined from the boundary condition. Although 
these equations have been integrated for large distances r, the boundary condition can 
be formulated only for small r — that is, at the body itself. However, here the bound- 
ary condition will be matched with the nonlinear solution at some distance R. If linear 
slender-body theory is used, the boundary condition means another approximation. But 
slender -body theory has here the advantage of short, compact solutions. To handle the 
problem in such a way means that the disturbances sent out by the body are assumed to 
run along straight characteristics (zeroth-order theory). At the distance R, the straight 
characteristics are matched with the characteristics resulting from the exact theory that 
are valid for large distances r. At first glance it seems to be a problem to find an 
appropriate matching distance R. But numerical testing showed that this parameter can 
be chosen in a broad range of values without essentially changing the numerical results 
(ref. 3). 

It is convenient to calculate the functions of integration by considering the angle of 
inclination 3 of the streamline. The inclination can be expressed by 


tan 3 = — - — (41) 

Uoo + u 

where v is a vertical velocity component and u^ + u is a horizontal velocity compo- 
nent. For small disturbances this equation can be approximated by 



M /3 9r p 9t7 


(42) 


where fi = \Jm 2 - 1 and T) = fir. 

If equation (42) is developed in perturbation and Fourier series, it can be rewritten 
as 


r* 

— 

r 

r- 


r - 

"I 

+ 3^2 cos 4 / + 

4^ + 4 2) J 2 c ° s w = < 

y 1 ) + y 2 > + 

y« + y*> 

2 cos i// + j 

<y> + y 2 ) 

2 cos 24 X 


, 

l 

L 



J 


(43) 
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By comparing equation (43) with equations (31) and (37), the functions of integration can 
be expressed in terms of cp^, which itself will be determined by slender-body theory. 
This will now be demonstrated for a lifting body of revolution. . . 


Lifting Body of Revolution 


The quantities <p^ can be taken from reference (3): 
^m 
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Qil 

Here A' = — ; a denotes the angle of attack, which is measured counterclockwise: and 
d£ . 

77 = /3R, where R is the distance from the body where the boundary condition will be 
satisfied. 

Thus the functions of integration are determined as 
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Within the limits of a small error, equations (45) can be integrated for a(z) = a(R). 

The integration functions G^ and G' 2 ^ are compared in figure 4 for a fighter 
airplane. The cross-sectional area distribution of that airplane is given in figure 5 for 
M = 1.6 Mach cuts. The negative areas are due to the fact that the engine stream tube 
exit areas are smaller than the intake. It can be seen from figure 4 that the integration 
functions satisfy the assumption of equation (16), since G^ » G® is the fundamental 
condition for the gas properties of first order, such as and w^, to be greater 

than those of second order. 


Pressure Signature 

The pressure coefficient is approximated by 


P - P 

A CO 

Ip u 2 

2 CO oo 


-2 — = -2 cos 9-1 
8x w<°> 


22 


The pressure Ap at the altitude z + r cos \p has been nondimensionalized with respect 
to the dynamic, pressure at the ground. By using the equation for the speed of sound at the 


the ground, 



o P 

— 7 T—> equation (46) can be rewritten as 
P gr 


Ap _ 

P gr 


n + 2 w<°> 

w 

n a 2 
gr 

w<°> 


cos 5 


1 


(47) 


To use this equation, the solutions for w and 3 have to be introduced. 


DISCUSSION 


Some pressure signatures will be given for a fighter airplane. The point of inter- 
section of the Mach line of the undisturbed flow through the nose of the body and a line 
parallel to the abscissa x at the distance r away from the body gives the origin of the 
plots of the pressure signatures, as shown in figure 6. The points of multiple values on 
the pressure signature indicate the occurrance of shocks. In figure 7 the pressure sig- 
nature is shown for the Mach number M = 1.5, the altitude of flight z = 10 000 m, and 
the speed of sound at the ground a gr = 300 m/sec at the distance r = 100 m away from 
the body. Here the factor a cos if/ has been varied. Naturally the pressure rise is 
greater at the bottom of the body (cos if/ = -1) than at the top (cos if; = 1) for an angle of 
attack a = -0.1. In figure 8 the decreasing pressure with increasing distance from the 
craft is demonstrated. In both figures only the bow shock has been calculated. A com- 
parison with the theory given in reference 2 is shown in figure 9. For this calculation a 
fairly good agreement has been obtained. 

At the caustic the speed of sound is of the same order of magnitude as the velocity; 
that is, w 2 - a 2 (r) « 1. For the cutoff point, this can be written as w 2 = a 2 (r) = s? 

- E(z + r cos if/). The locus of the cutoff can be calculated for given values of if/, z, a, 
and M(z): 


1 

a gr " m2 ( z ) &2 ( z ) 

I* — 

cutoff cos ^ 

E 


Because of the parabolic character of the flow, the pressure rise cannot be investigated 
at the cutoff point. At this point the assumed series, equation (16), breaks down. As can 
be seen from equations (31) and (32), the first-order angle of inclination decreases 

continuously and becomes zero at the cutoff point, whereas the first-order velocity com- 
ponent decreases at first and then increases, reaching a high value near the cutoff point. 
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Thus, at the cutoff point the first-order quantity is inconsistent with the assump- 
tion of equation (26), whereas the result for is compatible with this assumption. 

However, in the second-order calculation even the angle of inclination 3^, equation (37), 
contains singularities at the cutoff. However, so long as the Mach number is slightly 
supersonic, the pressure can be calculated. The pressure rise is shown in figure 10. 
First, a pressure decrease with increasing distance from the body can be observed. 

Then, ahead of the cutoff point, the pressure increases because of the focusing- and then 
augments very rapidly near the cutoff point. It should be pointed out that the most rapid 
increase of the pressure occurs within a distance of a few centimeters. For complete- 
ness the cutoff distance, at which the calculation breaks down, has been included in the 
figure. For this point, geometric acoustics give an infinite pressure rise. It is the . 
author's opinion that so long as the initial equations do not contain any damping mecha- 
nism, such as viscosity or entropy increase, the pressure rise at the cutoff cannot be 
expected to be finite. 

There are some uncertainties about how the pressure wave is reflected at the 
ground and at the cutoff point. Figure 11 (from ref. 5) shows the reflection of the pres- 
sure wave at the ground as measured from various tower microphones. It seems to indi- 
cate that the pressure wave is reflected in the same way as a wave at a solid wall; that is, 
the compression wave is reflected as a compression wave and the rarefaction wave as a 
rarefaction wave. The pressure signatures of the wave traveling toward the ground and 
the reflected wave can be superimposed. The abscissa of the reflected wave is x re ^ = -x, 
whereas the pressure signatures of the reflected wave are given by equation (47). How- 
ever, in this case no obstacles at the ground are taken into account; the ground is con- 
sidered to be plane and smooth. Figure 12 indicates the complexity of the boundary con- 
dition in the mixed-flow region of the caustic. It is not clear whether the atmosphere can 
be considered to behave as a "soft" solid wall at the cutoff point. Figure 13 (from ref. 5) 
shows some measurements at the cutoff point. Here the phenomenon of the U-shaped 
wave can be observed. The incident N-wave is changed at the focus with a component 
phase shift of -n/2 radians, whereas the reflected wave, after passage through the focus, 
again has an N-waveform. In first-order calculations the U-shaped waveform has not 
been obtained by the present theory. It is hoped to get an improvement by including the 
second order in the calculations. 

Although some progress has been made by Seebass (ref. 6) in calculating the sonic 
line at the caustic by using Tricomi's equation, this problem needs more investigation for 
the understanding of the complete behavior at the caustic and the cutoff point. 
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CONCLUSIONS 


The exact nonlinear system of partial differential equations for supersonic flow is 
solved for large distances in a nonhomogeneous atmosphere by using a perturbation 
method. The unknown functions were expanded in a perturbation and a Fourier, series. 
The system of differential equations was derived and solved for each order of magnitude 
and Fourier order. The integration introduced a function of integration which could be 
used to satisfy the boundary condition. To calculate the far field, the boundary condition 
was not satisfied at the body itself, but at some distance R away from the body. The 
matching point does not influence the numerical solutions in a broad interval of values 
of R. The example calculations give reasonable sonic-boom signatures. They are in 
good agreement with those obtained by previous theories. .s 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , February 14, 1973. 
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Wave normal cone 





50 m 



Influence of a cos \p on the pressure Figure 8.- The pressure signature 

re. M = 1.5; z = 10 000 m; R = 20 m; a cos \j/ = 0; M = 1.5; z = V 








Figure 9.- Comparison of present theory with theory of reference 2. Parabolic body 
shape, A = 0.047 t(x 2 - 2x^ + x^); a cos \p = 0; M = 1.5; R = 20 m. 
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Figure 11.- Measured sonic-boom signatures at various heights above the 
ground for an F-104 aircraft in steady, level flight at a Mach number 
above cutoff (M = 1.3) and an altitude of 10.26 km. (From ref. 5.) 







Figure 13.- Measured sonic-boom signatures at various heights above the 
ground for an F-104 aircraft in steady, level flight at cutoff Mach num- 
ber (M = 1.095) and an altitude of 10.26 km. (From ref. 5.) 
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